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(-H ■ Abstract 

We present a method to calculate nuclear matter properties in the superfluid phase. The method 
is based on the use of self-consistent off-shell nucleon propagators in the T-matrix equation. Such 
\& ' a complete treatment of the spectral function, is required below and around T c due to a pseudogap 

formation in the spectral function. In the superfluid phase we introduce the anomalous self-energy 
in the fermion propagators and in the T-matrix equation, consistently with the strong coupling BCS 
equations. The equations for the nucleon spectral function include both a contribution of condensed and 
scattering pairs. The method is illustrated by numerical calculations. Above T c pseudogap formation 
is visible in the spectral function and below T c a superfluid gap also appears. 
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1 Introduction 

2 . 

Nuclear matter calculations in the Brueckner approximation p], 0] do not consider a possible pairing in 
the ground state of nuclear matter ||, with the exception of Refs. Q. This approach is justified by the 
fact that the superfluid energy gap is expected to be small in normal density nuclear matter ||. When 
addressing the spectral properties of nucleons in medium, calculations using T-matrix approximation where 
performed p, pi M, p|. The T-matrix approximation for the two-particle Green's functions leads to a 
much stronger pairing than the Brueckner approximation. The T-matrix calculation and the BCS theory 
give the same prediction for the critical temperature 0. Brueckner type calculations of the effective 
in-medium nuclear interactions, permit to estimate approximate nuclear matter ground state energy , in- 
medium nucleon masses, and effective cross sections from the free nucleon-nucleon interaction. Although, 
discrepancies between different calculations and different parameterization of the nuclear interaction remain, 
a tremendous improvement over a simple Hartree-Fock approach is achieved. It is an interesting question, 
how to include the pairing force into this picture. 

The generalization of the T-matrix nuclear matter calculations below the superfluid transition temper- 
ature (T c ), requires by itself a procedure combining the ladder summation of the self-energy diagrams and 
a possible fermion pairing. On the other hand, the T-matrix scheme seems to be a natural starting point 
for the inclusion of superfluidity in nuclear matter calculations |9|, [10|. In this work we show that such 
a procedure is indeed possible, merging a T-matrix calculation of the normal fermion self-energy with a 
mean-field, frequency independent, superfluid energy gap. The interaction potential, which enters in the 
superfluid gap equation is the free nucleon-nucleon interaction. However, the scattering of nucleons and 
the pseudogap formation above the critical temperature reduces the value of the energy gap and of the crit- 
ical temperature, as compared to the BCS result. This can be achieved only using self-consistent off-shell 
propagators in the T-matrix diagrams [QT] . In such a way the formation of a pseudogap influences strongly 
the result for the T-matrix in attractive channels, shifting the appearance of a singularity in the T-matrix 
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at the Fermi energy to much lower temperatures. A consistent treatment of the superfluid transition re- 
quires then both the use of the off-shell nucleon propagators (below and above T c ) and of the anomalous 
self-energy (below T c ). To illustrate the formation of the pseudogap above T c , we present first the results of 
two calculations above T c (Sec. pf), a self-consistent T- matrix resummation, using off-shell propagators and 
a T-matrix calculation using on-shell nucleon propagators with Hartree-Fock single-particle energies. The 
results of that section indicate the necessity of using off-shell propagators in the T-matrix equation around 
T c) rather than a quasi-particle approximation. In section we derive the equation for the T-matrix self- 
energy and the gap equations, in the superfluid phase. A numerical example of a self-consistent calculation 
of the superfluid energy gap and of the T-matrix self-energy is presented is section |[ 

2 T-matrix equation in the Fermi liquid phase 

The nuclear matter can be described as an infinite system of interacting fcrmions. We consider a system 
of neutrons and protons interacting via two-body forces 
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H = Y,jd s x*Ux)(-^-n)*«(x) 



+ E l[d 3 x[d 3 y*l(x)*l,(y)V a ,^^(x,y)* (y)* a (x), (1) 

a' ,P' ,a,P 

Note that we define the energies with respect to the Fermi energy (chemical potential) \i\ also for real-time 
calculations. This will be helpful when dealing with the anomalous propagator in the superfluid phase. 
In the following we shall use equivalently the term zero energy and the Fermi energy. Nuclear forces can 
decomposed in (JST) (total angular-momentum, spin, isospin) channels [|13| 

^w(k, P )= E i l ~ 1 ' y^\k)v^\k^)y^' ipf . (2) 

(JST)U' 

The sum over the partial waves 1,1 is restricted to I = I = J for uncoupled states and to I, I = J ± 1 for 
the scattering of triplet states. 

One body observables of a systems of fcrmions at equilibrium can be described by the spectral function 



A(p,uj) = -21mG+(p,u>) 

-2ImS+(p,w) 



(u - p 2 /2m - RcS+(p, w) + m) 2 + ImE+(p, w) 5 



(3) 



where G denote the retarded(advanced) Green's function. We restrict the discussion to symmetric nuclear 
matter so that the Green's function is diagonal in spin-isospin indices. The T-matrix approximation 
consist of summing particle-particle and hole-hole ladder diagrams for the self-energy jlj, [lj, |lq| . The 
self-energy for in- medium propagators describes the dressing of the quasi-particles due to scatterings with 
other particles. In general one cannot approximate the spectral function by a S- function (quasi-particle 
peak) at some quasi-particle energy (& . The particles acquire a width, which enters self-consistently in the 
equation for the T-matrix. Moreover, more complicated structures in the spectral function seem to appear 
when the temperature is in a range of few MeV above the superfluid transition [pi . Using a quasiparticle 
ansatz for the spectral-function the authors of Ref. [p[ have shown that the spectral function has two (or 
sometimes three) peaks. This effect by itself invalidates the quasi-particle approximation for the T-matrix. 
On the other hand, it is possible to calculate the T-matrix and the corresponding spectral-function self- 
consistently using off-shell propagators in the ladder diagrams for the self-energy [O. The self-consistent 
calculation was performed in a simple model with separable interactions in the 5*-wave. In the present 
work we use this simple interaction for numerical examples. The self-consistent solution encompasses in 
a consistent way the non-Fermi liquid behavior of the spectral functions, and, thus, gives quantitatively 
different results for the critical temperature and for the temperature where the pseudogap appears, than 
the quasi-particle approximation. 

Off-shell propagators were used previously in a nuclear matter calculation by Jong and Lenske 112] . The 
spectral function in that work was approximated as a sum of a quasiparticle peak and a continuum part 
on the other side of the Fermi energy, the contribution to the two particle propagator with two continuum 



spectral functions was neglected to simplify the numerics. This neglected part could be calculated using 
the methods developed in the present work. On the other hand a mixed ansatz using a quasi-particle and 
a continuum part for the spectral function, like the one used in [Q or similar, is unavoidable at small 
temperatures. The authors of Ref [Q use a self-consistent scheme a zero temperature, but do not consider 
pairing. However, even if the expected pairing is small it is more natural to use a scheme allowing for the 
creation of pair condensate, since the normal Fermi liquid is unstable. 

The T-matrix for a system with a two-body interaction V(p, p ) is defined as GL Qq, Eq, ftjl : 

d 3 k f d 3 c 
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< k|^(P, W )|q X q|T* a/5 (P,uO|p > , (4) 

where the disconnected two-particle propagator is : 



< P |e±(P, W )|p >= (2^) 3 <5 3 ( P -p)J^-J. 



du) f duj 

G < (P/2 + p,oj" -u>')G < (P/2-p,u}') 
-G > (P/2 + p,co" -w')G > (JP/2-p,w'))/(w-w" ±ie) , (5) 



with 



G<(p,u) = iA(p t u)f(u) 

G>( P ,u;) = -M(p,a;)(l -/(«)) , (6) 

and 

is the Fermi distribution. The in-medium T-matrix equation can be decomposed in partial waves if an 
angle averaged two-particle propagator (G)(i is used : 

(< p|g ± (P, w )| P >) n - J ^ < p|0 ± (P,«)|p' > , (8) 

where we average over the angle between P and p. 

<k\g ± (P,Lo)\k>< k\T^ T) ± (P,uj)\p > , (9) 

In the course of this work we shall discuss two different approximation schemes for the intermediate two- 
particle propagator in the T-matrix calculation. The formulas for the in medium T-matrix approximation 
in nuclear matter with partial wave decomposition can be found in Ref. |17| . 

• Self-consistent solution 

The self-consistent solution uses off-shell nucleon propagators in the ladder diagrams (Eq. Q). The 
imaginary part of the self-energy is defined by the T-matrix in the following way 

(JST)l J J y ' 

< (p - k)/2|ImT// ST) + (|p + k|, uj + w')[(p - k)/2 > 



A(k,u>)(f(u')+b(u + u'j) , (10) 



where 

& M = ^w~ x ( n ) 

is the Bose distribution. The real part of the self-energy is the sum 

ReE(p,c J ) = S ffF (p) + S d (p,w) (12) 

of the Hartree-Fock term 

d 3 k 



w-iEP^ii/t/H 



(JST) 



V^ JST) (|p - k|/2, |p - k|/2) A(fc, w)/(«) , (13) 

and the dispersive contribution to real part of the self-energy 

„ , . _ f du> — ImE + (p, u ) , , 



The solution of Eqs. (H), (10), rtlq), (|l4|), and (0) is obtained by iteration with the constraint 

£/|^ (P,W)/(U ' ) = P ' (15) 

where p is the assumed density of the nuclear matter. Numerical methods used in the solution of the 
self-consistent T-matrix equations are presented in an appendix. 

• Quasi-particle approximation 

The quasi-particle solution is obtained by using mean-field propagators in the T-matrix equation. 

G ™f^ = ^±i-e> (16) 

where £ p = u p — fi is the mean-field energy measured with respect to the Fermi energy, and the mean- 
field quasi-particle energy is of course lv p = p 2 /(2m) — T,hf{p)- More generally, the self-consistent 



P = A I 7^3 /&>)• ( 17 ) 



real part of the self-energy (12) can also be used to define the quasi-particle energy. The density is 
given by 

d 3 p 

Analogously the expressions for the self-energies in the quasi-particle approximation take the form 

ShfM , E(2r+1)(2J+1) | 

(JST)l y ' 

^ (JST) ((| P -k|/2,|p-k|/2)/(6), (18) 

d 3 k 



Im£+(p,^) = ^- Y, (2T + 1)(2J +!)/"■ 



■ (2tt) 

(JST)l " v ; 



< (p - k)/2|ImT// ST) +(|p + k|, w + &)|(p - k)/2 > 

(/(&) + &(& + <"))• (19) 



The dispersive contribution to the real part of the self-energy can be obtained from the dispersion 
relation ( |l4| ) . Note that there is no need for an iteration procedure in the solution of the T-matrix 
equation in the quasi-particle approximation, when using Hartree-Fock single particle energies, as in 
this work. Once the mean-field self-energy is calculated, the T-matrix can be obtained as 

.-.-.i / S' I i I _ , ' - .; / S' D ' ^~ • > / I '"' 

* ' I I •■ ' ' '- 11 " - - ' 

l" 

(1 - f(U f(Uh < ^JST) ± (p>w) | p ' > ; (20) 



Pir^ Hp^)\p >=^ T W')+£/|^rW) 



Pl , 2 = P/2±k. 



For numerical convenience, we calculate the angle average separately in the numerator and in the 
denominator of the two-particle propagator in the above equation. The spectral function is then 
calculated using Eq. ([l9|), the dispersion relation (fL4|), and the definition (||). If the spectral function 
presents a quasi-particle peak around the energy £ p , the quasi-particle approximation can be trusted. 

Of course an iterative improvement of the single particle energies using the real part of the self-energy 
obtained in the T-matrix or G- matrix approximation is possible M, pj. Such an iterative scheme 
requires in principle the introduction of BCS-like single particle energies in the presence of pairing. 
This can be calculated in a quasi-particle version of the formalism used in Sect, m However the actual 
implementation of this approach is quite tedious and we restrict ourselves to the Hartree-Fock single 
particle energies and to the normal Fermi-liquid phase for the the quasi-particle approximation. 

A common property of all the T-matrix approaches to the nuclear matter is the appearance of strong 
pairing in channels with attractive interactions. As the temperature is lowered a singularity builds up in 
the T-matrix at zero total momentum and energy. For the quasi-particle scheme it leads to the Thouless 
criterion for T c , which is equivalent to weak coupling BCS |J. The inclusion of self-energy corrections in 
the normal fermion propagators generalizes this criterion, giving more general equations for the critical 
temperature. In section [| we shall derive equations valid also below T c . 

3 Pseudogap formation in finite temperature nuclear matter 

Numerical results presented in this, and subsequent sections were obtained using a Yamaguchi separable 
potential of rank one |lq ]. The momenta of all the nucleons are restricted to \p\ < 700 MeV. The coupling 
strengths of the original Yamaguchi parameterization are slightly adjusted to reproduce the scattering 
length in vacuum as obtain from the Yamaguchi potential without cutoff in momenta of nucleons 

V(k,p) = Yl ^g(k)g{ P ) , 

Q = 1 S , 3 Sl 

9(p) - pnb* ' (21) 

with Ai So = -.841 GeV 2 , Aa Sl = -.8591 GeV 2 , and /3 = 285.9 MeV . For large momenta the effect of the 
cutoff is significant. Nevertheless, the following calculation is a helpful illustration, presenting all the basic 
phenomena related to the superfluidity in nuclear matter, and describing qualitatively the behavior around 
the Fermi energy. The binding energy of the deuteron and the critical temperature (in the quasi-particle 
approximation) are changed within few percent by the introduction of the cutoff. 

All the results are calculated for nuclear density p — A5p . At higher densities, a resonance appears 
in the two-particle Green's function below the Fermi energy. This happens for temperatures above T c , 
according to the Thouless criterion. This statement is strongly dependent on the assumed interaction and 
the resulting deuteron wave function. This phenomenon is not related to the main subject of this work. We 
postpone the discussion of this topic, using more realistic parameterizations of the momentum distribution 
of the deuteron, to another work and restrict ourselves to low density nuclear matter in the following. The 
Yamaguchi potential has also the practical advantage, that the critical temperature and the energy gap 
are relatively (somewhat unrealistically) high [g|. The strength of the pairing simplifies the numerics, and 
is helpful in illustrating qualitatively the physical effects. In Ref. [ll]] we have studied the self-consistent 
T-matrix approximation at normal nuclear density and high temperature. However we were not able to 
perform similar calculations at small temperatures (around and below the critical temperature) with the 
limited resolution in energy, because the single-particle width becomes very small at the Fermi energy (See 
appendix). 

3.1 Quasi-particle approximation 

The spectral function in the quasi-particle approximation has been calculated in several works |5J, |3|, |l9|, |20j , 
in the T-matrix or Brueckner schemes. Calculations using the T-matrix approximation for the fermion self- 
energy show very strong pairing in channels with attractive interactions. At a critical temperature a pole 
appears in the T-matrix. For the interaction here chosen, it appears first in the in the 3 5i channel. It signals 
the formation of neutron-proton bound states (with zero binding energy, since the pole in the T-matrix 
appears at zero total momentum and at the Fermi energy) |j, g , for even lower temperature it is preferable 
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Figure 1: Imaginary part of the self-energy 
— ImS(p, tu) as function of the energy for p — and 
p = 175 MeV (solid and dashed lines respectively), 
calculated in the quasi-particle approximation for 
the T- matrix at p = .54p , and T = 5.1 MeV 
(slightly above T c =5.03). 



Figure 2: Spectral function A(p,u>) for several val- 
ues of momentum as function of energy for the same 
parameters as in Fig. pi The lines correspond to 
p = 0, 70, 140, 210 and 350 MeV, according to the 
position, from left to right, of the largest peak in 
the spectral function. 



for the system to condense into nucleon-nucleon pairs instead of remaining in the Fermi-liquid phase. At 
temperatures slightly above T c , the T-matrix is strongly peaked for small values of total momentum and 
small energy. It is a precursor of the singularity in the T-matrix at the critical temperature. One of the 
consequences of this singularity is a very unusual form of the imaginary part of the self-energy (Fig. PP) . The 
imaginary part of the self-energy Im£ + (p = 0, ui) has a peak for u — —£ P =o- Not surprisingly, the spectral 
function exhibits also a nontrivial structure, with two or more peaks (Fig. 0). The spectral function is not 
of a Fermi-liquid type. 

The appearance of two peaks in the spectral function above T c was broadly discussed in high T c supercon- 
ductivity under the name of pseudogap formation pi, 22, 23L E4 GB, E6|. The quasi-particle approximation 
represents only the first iteration in the self-consistent calculation of the nucleon spectral function. Thus, 
it fails to describe the feedback of the pseudogap on the value of the T-matrix in the vicinity of the Fermi 
energy. This leads to quantitative difference between the self-consistent and the quasi-particle calculations. 
The critical temperature is 1.63 MeV for the self-consistent calculation and 5.03 MeV in the quasi-particle 
approximation. The form of the spectral function is different in the two approximations Jy] . The strong dif- 
ferences between the quasi-particle approximation and the self-consistent solution indicate that for strongly 
interacting systems (having bound or almost bound states in vacuum) the usual approximation of the spec- 
tral function by a quasi-particle peak fails already above T c . Another difference already mentioned in |11 
is that the imaginary part of the self-energy has a singularity in the quasi-particle approximation; in the 
self-consistent calculation this singularity is smeared out. In Fig. [|we plot the position of the maxima of 
the peaks in the spectral function as function of the position of the quasi-particle pole C, p ; 



Cp 



2m 



ReS(p,C P ) + M = 



(22) 



Clearly two peaks are visible in the spectral function on two sides of the Fermi energy. The dominant peak 
follows approximately the position of the quasi-particle pole. However, a subdominant peak is found on the 
other side of the Fermi energy, recalling the effect of an energy gap in the superfluid phase. The effect was 
named a pseudogap, since the energy gap in the spectral function occurs without pairing condensate and 
the superfluid density is zero. Close to the Fermi energy the subdominant peak disappears. The position 
of the quasi-particle peak (solution of Eq. p2[ ) is relatively smooth as a function of momentum, in both the 
quasi-particle approximation and in the self-consistent solution. There is no sharp variation of the effective 
mass near the Fermi energy. 
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Figure 3: The positions of the largest peaks of the 
spectral function at positive and negative energies 
(dotted and dashed lines) as function of the energy 
of the quasi-particle pole, for the same parameters 
as in Fig. 111. The solid lines denote the two asymp- 
totic branches in the BCS solution E p = ±£ p and 
the Fermi energy E p = 0. 



Figure 4: Spectral function for several values of mo- 
mentum as function of energy. The lines correspond 
to p = 0, 70, 140, 210 and 350 MeV, according to the 
position, from left to right, of the largest peak in 
the spectral function. The results were obtained in 
the self-consistent scheme at T = 1.7 MeV (slightly 
above T c = 1.63 MeV) and p = .45p - 



3.2 Self-consistent solution 

The self-consistent solution has a similar behavior around the critical temperature as the quasi-particle 
approximation: at the critical temperature the T-matrix has a singularity at zero total momentum and 
energy, and for temperatures above T c , the T-matrix is strongly peaked for small values of energy and 
total momentum. This leads to a pseudogap structure in the spectral function (Fig. Eh. This pseudogap 
structure is again similar to the superfluid gap. Two peaks of the spectral function are located on both 
sides of the Fermi energy. The dominant peak (with the larger weight) is the one whose energy is the closest 
to the energy of the quasi-particle pole. Please note that for momenta close to the Fermi momentum there 
is only one peak in the spectral function (the curve corresponding to p = 210 MeV in Fig. Q). In Fig. g 
the dispersion of the two peaks of the spectral function is presented. The dispersion of the maximum of 
the two peaks is very reminiscent of the superfluid dispersion relations Eu — ±vCfc + A(fc) 2 , where A(fc) 
is the pseudogap. At the Fermi energy the separation between the two peaks should be equal to twice 
the value of the pseudogap. However, close to the Fermi energy only one, dominant peak of the spectral 
function survives. Its width is small close to the Fermi energy, and the position of its maximum is given 
by the energy of the quasi-particle pole £ p . 

It is difficult to devise a self-consistent approximate scheme, using a two-pole ansatz for the spectral 
function. Usually such a procedure would require, when iterating Eq. (IQ), the introduction of additional 
poles in the spectral function. A self-consistent approximate solution in the pseudogap region can be 
obtained only in the pairing approximation p3[ , i.e. using a mean- field propagator on the right-hand 
side of Eq. (|l0|). If additional assumptions are made for the T-matrix in the pseudogap region, a self- 
consistent approximation scheme with superfluid- like ansatz for the spectral function can be obtained p3[ . 
However, the self-consistent solution indicates (Fig. H), that the pseudogap picture |23[ of single-particle 
excitations above T c is not entirely correct. Moreover, the energy of the quasi-particle pole is different 
from the Hartree-Fock energy. The imaginary part of the self-energy and its energy dependence are usually 
non-negligible. 

In Fig. O the integrated weights of the two peaks in the spectral function are presented. The weight of 
the peak is obtained by integrating the spectral function on the same side of the Fermi energy. Although 
qualitatively the weight of the two parts of the spectral function follows the superfluid result iV,(fc) = 

| ( 1 ± j*- I , we could not fit the weight N^ using this law. This is not surprising since the two peaks 
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Figure 5: The positions of the two peaks of the 
spectral function (dotted and dashed lines) as func- 
tion of the energy of the quasi-particle pole, for the 
same parameters as in Fig. El The solid lines denote 
the two asymptotic branches in the BCS solution 
E p = ±£p and the Fermi energy E p = 0. 



Figure 6: The weight of the peaks of the spectral 
function at positive and negative energies, for the 
same parameters as in Fig. EJ 



are very broad, and, moreover close to the Fermi energy where the effect of the pseudogap should be the 
clearest, the two peaks merge into one. The sum of the weights of the two peaks gives of course 1. 

The nuclcon momentum distribution is presented in Fig. |?]. The solid line was obtained using the 
self-consistently calculated spectral function 



,/(/»= I^A(p, U )f(u>) 



(23) 



The dotted line denotes the result of the quasi-particle approximation for the spectral function 

n{p) = f(Cp) , (24) 

where ( p denotes the position of the quasi-particle pole, i.e the solution of Eq. (f22|). Nucleon scattering 
(nonzero imaginary part of the self-energy in A) leads to a considerable broadening of the momentum 
distribution. This is another indication that the quasi-particle approximation is not reliable for nuclear 
matter around the critical temperature. 



4 Pairing and the T-matrix approximation 

In the superfluid phase the nuclear matter develops a superfluid energy gap and non-zero anomalous 
fermion propagators. The usual BCS theory is limited to mean-field approximations to the self-energy 
and to the superfluid gap. Formally T-matrix equations can be written for the normal propagators and 
anomalous propagators ||27| , [22| . As before, T-matrix enters in the calculation of the self-energy in place of 
the interaction potential. But also a generalized T-matrix, between the anomalous propagators enters in 
the calculation of the anomalous self-energy. 

The resulting equations are, however, very difficult and no solution in this scheme has been presented. 
In the following we will present the simplest generalization of the T-matrix and BCS schemes, allowing for 
the calculation of the spectral function in the superfluid phase. It is known that there is a link between the 
T-matrix and the BCS theory at the critical temperature. The critical temperature in the BCS theory is 
related to the appearance of a singularity in the T-matrix at zero energy H . It is a manifestation of more 
general relation between the T-matrix fermion scattering and pairing, and the BCS gap equation. This 
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Figure 7: Nucleon momentum distribution obtained 



from the spectral function (Eq. 23 ) (solid line) and 
using the quasi-particle pole approximation for the 
spectral function (dotted line), for the same param- 
eters as in Fig. 0. 



Figure 8: Diagrammatic representation of the T- 
matrix equation in the superfluid phase. The thick 
line denotes the off-shell normal fermion propa- 
gator, the thick line with a dot denotes the full 
fermion propagator including the anomalous self- 
energy, the wavy line denotes the interaction and 
the box stands for the T-matrix. 



relation was observed long time ago by Kadanoff and Martin jlC| ]. They proposed to use one mean-field 
fermion propagator and one BCS two-pole fermion propagator in the ladder of the T-matrix diagramR : 

d 3 k f dtu f duj . , 

j^rJ^J -^rw 7 *(p>K) 

A(P/2 + k, J - uj")A s (P/2 - k, w")(l - f(u - u") - f(u>")) 



lu — to' ± it 



<k|T± Q/J (P,a,)|p'>, 



(25) 



where A is the spectral function without anomalous self-energy, which we will call in the following normal 
fermion spectral function, and which takes the form 



A(p,w) = 27nJ(w-&) , 



(26) 



in the quasi-particle approximation, and A s is the spectral function with anomalous self-energy, which we 
call full fermion spectral function in the rest of the paper (Fig 0) . In the mean-field approximation it takes 
the usual BCS form 



A.(p,w) = 2nl^±^S^ - E k ) + %-^5{u + E k ) 



2£ k 



2£ k 



where 



s k = ,/e k + A(k)2 



Using ( |25| ) the gap equation can be written as 

d 3 p 



Y,l^k< PlHtfJ^CP = 0,- = 0)|p > A Q/3 (p ) = , 

a/3 V ' 



(27) 
(28) 

(29) 



1 We use the real time formalism also in the superfluid phase J28|, with energies measured with respect to the Fermi energy. 



the relation between A and A a p is given later. The limit of the above equation when A = in the T-matrix 
equation, gives the Thouless criterion for T c . The meaning of Eq. ( p9| ) is, however, more general. It states 
that the gap equation is equivalent to requiring that the singularity of the T-matrix remains at zero energy 
and momentum for all temperatures below T c . The principle of conserving approximation, in the sense of 
Kadanoff and Baym |14| , requires the use of the normal fermion mean field propagator (without anomalous 
self-energy) on the right hand side of Eqs. ( |l0| ) and (|l3|). 

In the above equations we have assumed that the ground state is time- reversal invariant. Thus, the full 
fermion propagator G s is still diagonal in spin-isospin indices 129, BOl. However, the full spectral function 
and the full propagator depend in general on the direction of the momentum. Also the T-matrix would 
then depend on the direction of the total momentum of the pair. 

In principle the pairing approximation of Kadanoff and Martin can be used to calculate the spectral 
properties of the nuclear matter in the superfluid phase using the quasi-particle approximation. However, as 
we have seen, the quasiparticle approximation becomes unreliable in the case of strong attractive interaction 
between fermions, with bound or almost bound states in vacuum. As discussed in models of high T c 
superconductivity, the crossover region between the BCS pairing and the boson condensation, cannot be 
described using mean-field fermion propagators for the calculation of the T-matrix. The same is true in 
nuclear matter for with the Yamaguchi interaction in the T-matrix approximation, because of the pseudogap 
formation in the spectral function (Sec. g). At higher densities, in particular at the saturation density, 
the pairing gap is expected to be smaller and the BCS theory is correct. Certainly calculations using more 
realistic interactions for several densities, including the modification of single-particle energies and finite 
quasi-particle width, would clarify at which densities, if any, the crossover-like behavior between the BCS 
pairing and the boson condensation occurs. One way to generalize the Kadanoff-Martin approach was 
proposed by the Chicago group [^3[ |24j. From the observation that a pseudogap in the spectral function 
forms above T c , the authors of Refs. [g3| g4j proposed to use one mean-field and one full propagator 
in the T-matrix equation. The full propagator would include self-consistently the normal self-energy in 
the T-matrix approximation (defined by (nq) and (J19|)) , and, below T c , also the anomalous self-energy ( 
defined by the superfluid gap (plj)). In actual calculations, the authors of J2|, |24j used an approximate 
BCS-like form also the normal self-energy in the vicinity and below T c . The spectral function is then 
given in an analytical form by the superfluid gap and by a pseudogap, defined by the singularity of the 
T-matrix for small energy and small total momentum. This approach represents an improvement over the 
quasi-particle approximation, because it takes into account the pseudogap formation in the calculation of 
the T-matrix. Also its a very nice illustration of the mechanism of the pseudogap formation and of the 
two-fermion excitations close to the Fermi energy Q. The pseudogap, as well as the superfluid gap, causes 
the formation of an energy gap in the two-fermion excitations around the Fermi energy, very similar to 



the one postulated in W. However, the self-consistent spectral function calculated in Sec. 3.2 shows only 
partly a pseudogap. The position of the subdominant peak in the spectral function does not follow the 
relation expected from a pseudogap energy dispersion relation. Also close to the Fermi energy only one 
peak in the spectral function is present. One could use the full spectral function instead of the pseudogap 
approximation for one of the propagators in the T-matrix equation. Thus, the correct dispersion of the 
peaks in the spectral function would be taken into account self-consistently. However, it turns out that 
the numerical cost of a calculation in the pairing approximation with one off-shell propagator and of a full 
self-consistent calculation is similar. Also the use of the mean-field propagator in the calculation of the 
self-energy (|l9|), as required in the pairing approximation, leads to a singularity in the self-energy at T c . 
No such effect is present in the self-energy from the self-consistent calculation |11| . 

Therefore, we shall seek a generalization of the self-consistent T-matrix equations in the superfluid phase, 
including the off-shell fermion propagators in the T-matrix equation and in the self-energy calculation. Since 
the T-matrix equation in the pairing approximation is related to weak-coupling BCS equation [|10[ , we will 
use a mean-field approximation for the anomalous self-energy. Indeed the superfluid gap equation (Fig O) 

A Q/3 (P) = - E / £ / (03 *W ((P - k ) A (P - k)/2)B^, (k, W )/M , (30) 

a 

can be related to the T-matrix equation, using strong coupling BCS equations, with frequency independent 
gap parameter || || (Fig |l^) 

G+(p,c) = G+(p,uO - ~X>^(p)Fi+(p,uO (31) 

a/3 

Fj[/(p,u;) = G-(-p,-«)G+(p,w)At / ,(p) , (32) 

10 



-<3> 



(a) 



"«— »- 



(a) = cL~$ 



4 ► 



(b) 



^ 




■*<£>*- 



(c) 



Figure 9: Diagrammatic representation of the gap 
equation (p3CJ) (first line). The thick line with two 
arrows denotes the anomalous propagator and the 
circle with the letter A denotes the pairing gap, 
other symbols as in Fig. H. The equality in the 
second line, is equivalent to Eqs. d2J3), (pi). 



Figure 10: Diagrammatic representation of the BCS 
Dyson equation, (a) represents the resummation of 
the anomalous self-energy in the nucleon propaga- 
tor in the supcrfiuid phase (Eq. |3l]), (b) represents 
the equation for the anomalous propagator (Eq. p2|) 
and (c) represents the usual resummation of the 
self-energy (the circle with the letter £). Symbols 
as in Figs. @, §. 



where G s is the full fermion Green's function, including the anomalous self-energy (G is reserved, from 
now on, to the normal Green's function including only the normal self-energy). The anomalous propagator 
spectral function is 

B aP (p,u) = -21xxiF+ p (p,oj) . (33) 

Using (f32|) in the gap equation (p0|), and using the spectral representation of the propagators G and G s 
one obtains 

, . f d 3 k f dui f du> rr , , . 

A(P/2 + k,a; - cu')A s (P/2 - k,u)(l - f(oJ - J) - f(u')) 

w 

A 75 (k)=0. (34) 

Thus, the gap equation for the superfluid gap is equivalent to the condition that A is the zero eigenvalue 
vector of the real part of the inverse T-matrix (Eq. 03), i.e Eqs. (p9|), (B4) and ( p0[ ) are equivalent. However, 
here the T-matrix is given by the expression (pa), with the spectral functions A and A s being obtained 
from the solution of @ (A(j),u>) = -2TmG + (p,u>), A s (p,uj) = -2ImG+(pw)). 

The set of equations (p5|), (pl|), and (pp|), defines a consistent approximation scheme for the nucleon 
self-energy below T c . In this approach the singularity of the T-matrix is always at zero energy and zero 
total momentum, the same is true in the BCS approach p0|. However here we include the contribution of 
the scattering and resonant pairs to the self-energy, through Eq. ([H]), as well as of the condensed pairs, 
through the anomalous self-energy. In the region of intermediate coupling strength, between the weak 
BCS and the Bosc condensation, the contribution of the T-matrix to the normal self-energy is important 
p7|, pi| , |23| , gjj , and leads to a pseudogap formation. 

There is one more technical obstacle before using this scheme in actual calculations. In the superfluid 
phase the full propagator G s (p, u>) depends on the direction of the momentum p, and the T-matrix depends 
on the direction of the total momentum. This dependence introduces tremendous difficulties in the nu- 
merical solution of the T-matrix ( poj ) and the BCS Dyson equation (|3l|). Also in the T-matrix calculation 
we have used partial-wave expansion, enabled by the angle averaging of the two-fermion propagator in 
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the ladder diagrams. On the other hand, the solution of the superfluid gap equation, which is equivalent 
to the calculation T-matrix equation at zero total momentum, requires no angular averaging of the two 



propagators in (32). In the BCS theory an approximate decoupling of different angular momentum gaps 
is obtained p3|. In that case the gap equation in each partial wave approximately decouple, i.e. we can 
use a different propagator G s in each of the partial wave gap equations. In nuclear physics the situation 
is more complicated due to the coupled partial waves in the T-matrix. If we include the coupled partial 
waves, we cannot use the partial wave decoupling of the gap equations (e.g. 3 Si — 3 Z?i). Also in the nuclear 
matter we can in general have two superfluid gaps for the same partial wave but with different total angular 
momentum and isospin (e.g. So and Si). Moreover, we want to use the angular averaged two-nucleon 
propagator in the T-matrix equation. Thus, the corresponding gap equations in different channels will 
be coupled through the dependence of the full fermion propagator G s on the anomalous self-energy. In 
practice only the gap corresponding to one channel would be non-zero, the one with the strongest pairing 
at a given temperature and density. The T-matrix in other channels would have the singularity pushed 
away from the Fermi energy due to the energy gap originating form the paired channel. Of course this 
does not exclude phase transitions between pairing in different channels, when changing the temperature 
or density. 

In order to simplify the description of the full fermion propagator we shall use an angle averaged 
anomalous self-energy in (|3l]). This approximation simplifies also the gap equation |30|], in particular the 
superfluid gap for non-zero angular momentum can be taken as real and independent of the projection of 
the angular momentum [M . Of course such an assumption for the superfluid gap is only an approximation, 
and in some cases the true ground state of the nuclear matter may be described by more complicated angle 
dependent anomalous propagators. The example of the liquid Helium teaches us that the structure of 
the anisotropic superfluid order parameter can be very rich |pq| . In particular the real order parameter 
would not lead to the state of lowest energy. However, it is sufficient for consistent and stable T-matrix 
calculation of the nuclear matter. It would be extremely difficult to implement the most general angle 
dependent superfluid gap and full fermion propagator in the T-matrix calculation. 

As a result the self-consistent T-matrix can be decomposed in partial waves 

<P\T^ ST) ± (P^)\p >= V^ ST \p,p) 
,A(|P/2 + p|, «' - u")A s {\T>/2 - p|,w")(l - /(u/ - w") - /(«")) 



u> — lu' ± it 



The full fermion Green's function can be written as 



n 



,(JST) ±, 



<k\T^ I}± (P,co)\p'> . (35) 



Gf(p,u>) = 7^-7 , 1 : A2 , ,„ , r , (36) 

G+(p,u) 1 + A 2 (p)G (p,-w) 

where the normal fermion propagator is given by the standard expression and A(p) is the angle averaged 
total energy gap. The self-energy E is given by Eqs. J13] ) and (|10|), with the normal spectral function A 
on the right hand side. The set of gap equations in (JST) channels takes the form 

2, 



^'-?/5/^/5?Tfc'> 



A(k,uj- uj')A s (k,u')(l - f(uj - J) - f(uj')) A (JgT) ( - fc s , 37) 

U) l' ' 

The existence of a non-zero solution of the above equation in a particular channel is equivalent to the 
presence of a singularity in the T-matrix in that channel at zero total momentum and zero energy. The 
total, angle averaged energy gap is given by an incoherent sum of contributions in different channels 



A2 ^ = T- E (2T+l)(2J+l)A\ JST \pf 



«7T 



(JST)l 



- -±-(2T+l)(2J+l)J2^ JSTh (p) 2 , (38) 



87T 
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where in the last equality we have used the fact that only one channel (JST) S gives nonzero gap. It is the 
channel which gives the larger value of A(p). One can notice that using mean-field spectral function in 
(|37|), one recovers the BCS gap equations in (JST) channels from Q. The final equation of the calculation 
scheme is the expression for the nuclear density 

"= 4 /^/w fi4 * (p,w)/(fa,) - (39) 

It is instructive to write explicitly the full fermion spectral function in terms of the energy gap and the 
self-energy 

A s { P ,lu) = -2((w + £ p +Re£ + (p,-w)) 2 Im£ + (p,w) 

+Im£+0, -uj)A 2 (p) + (Im£+(p, -w)) 2 ImE(p,w) ]/ 

( ((w - & - Re£+(p, w)) (w + £ P + ReS+( P , -w)) 

\ 2 
-ImE+(p, w)ImS + (p, -w) - A 2 (p)J 

+ (lm£ + (p, w)(w + C P + ReS+(p, -w)) 
+Im£+(p,-u;)(a;-£ p -Re£ + (p,w))) 2 j . (40) 

It is clear that in order to find the spectral function for energy w, one has to know the self energy for 
energy u) and — u>. Therefore, the numerical solution must be performed in a symmetrical interval around 
the Fermi energy (a; = 0). 

5 Nuclear matter in the superfluid phase 

We solve numerically the set of equation ( |35|) , (|Tc|) , (pOJ) , (p^), with the constraint (^) for the same 
interaction as used in Sec. 0. Eqs. ( p9| ) and ( p7[ ) are solved simultaneously for the superfluid gap A(fc) and 
for the chemical potential //. Below T c the gap equation has a nontrivial solution, corresponding to the 
superfluid phase. For the case of the Yamaguchi interaction, the pairing occurs always in the 3 S\ channel. 
The T-matrix has a singularity at zero momentum and zero energy in this channel, for all temperatures 
below T c . The T-matrix in the 1 Sq channel has no singularity but is strongly peaked near the Fermi 
energy. For more realistic interaction we expect also pairing in the lg channel and pairing with higher 
angular momentum at larger densities. The calculation gives two spectral functions, the normal spectral 
function A{p,ui), without the anomalous self-energy, and the full spectral function A s (p, u>), obtained from 
the BCS Dyson equation (|3l|). We have observed that due to strong attractive interaction, and due to 
the appearance of the singularity in the T-matrix, a BCS-like two peak structure appears in the spectral 
function A in the pseudogap region. The same is true below T c (Figs. |lj and |IJ). The spectral function 
without anomalous self-energy is very similar as above T c (Fig. ^). Indeed the imaginary part of the 
self-energy (Fig. [ll]) is almost indistinguishable for the two temperatures, except near the Fermi energyg. 
In the full spectral function another mechanism is also responsible for the double peak structure. For all 
momenta the full spectral function A s has two peaks situated on both sides of the Fermi energy. Close 
to the Fermi energy, the superfluid gap in the spectral function is clearly visible (Fig O). Far from the 
Fermi energy the modifications of the dominant peak are negligible. The subdominant peak of the spectral 
function is only slightly modified by the anomalous self-energy. This is what is expected form the usual 
BCS theory, since the weight of the subdominant superfluid peak in the spectral function i ( 1 — 'j=P- ) is 



2 V Ek , 

becoming small far from the Fermi energy. 

The dispersion of the maximum of the peaks of the spectral functions A and A s illustrates the formation 
of the superfluid gap (Fig. 09) . On the abscissa we plot the energy of the quasi-particle pole in the normal 



2 As expected the scattering width at the Fermi energy decreases with the temperature. 
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Figure 11: The imaginary part of the self-energy, 
calculated in the self-consistent T-matrix approxi- 
mation as function of the energy, for the same pa- 
rameters as in Fig. 0. The solid and dashed lines 
denote the results at p = and p = 175 MeV re- 
spectively. 



Figure 12: The spectral function A(p, w) (solid 
line) and the spectral function with anomalous self- 
energy A s {p,ui) (dashed line) as function of energy 
at p = 140 MeV. The results were obtained at 
T = 1.4MeV, and p = A5po, corresponding to a 
superfluid gap of A(0) = 5.1 MeV. 
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Figure 13: Same as in Fig. nj but for p = 210 MeV. Figure 14: Same as in Fig. pj but for p = 280 MeV. 
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Figure 15: The positions of the two peaks of the 
normal spectral function (dotted and dotted lines) 
and of the full spectral function (dash-dotted lines) 
as function of the energy of the quasi-particle pole, 
for the same parameters as in Fig. [13. The solid 
lines represent the two asymptotic branches in the 
BCS solution E p — ±( p and the Fermi energy E p — 
0. 



Figure 16: The weight of the peaks of the normal 
spectral function at positive and negative energies 
(dashed and dotted lines) and of the peaks of the 
spectral function with anomalous self-energy (solid 



lines), for the same parameters as in Fig. 12 



fermion Green's function G. This energy coincides approximately with the maximum of the dominant 
peak in the spectral function A. This is especially well fulfilled in the vicinity of the Fermi energy. The 
subdominant peak in the spectral function A is visible for energies further away from the Fermi energy, 
also its energy does not follow the branch E p — ~( p . Like for the case above T c , the pseudogap in the 
spectral function is only qualitatively similar to a superfluid gap. Close to the Fermi energy, the spectral 
function A is similar as in the normal Fermi liquid. On the other hand, the positions of the peaks of the 
full spectral function A s follow the dispersion relation expected for a superfluid matter, especially close to 
the Fermi energy. An energy gap is clearly visible for C, p ~ 0. Away from the Fermi energy, the positions 
of the two peaks approach E p = £ p for the dominant peak and E p = — £ p for the subdominant peak in the 
full fermion spectral function. For \( p \ > 10 MeV, the positions of the peaks do not follow the BCS energy 
branches any more. However, one can observe that the dominant peaks in the normal fermion spectral 
function and in the full fermion spectral function are very close away from the Fermi energy, as expected. 

In Fig. nq are presented the weights of the peaks in the spectral functions A and A s as function of the 
energy of the quasi-particle pole of the spectral function A. The peaks in the two spectral functions have 
qualitatively similar behavior. However, the weight of the dominant (subdominant) peak is smaller (larger) 
in the full spectral function, than in the normal spectral function. It means that part of the weight of the 
spectral function is shifted from the dominant to the subdominant peak, as expected from the presence of 
an additional mechanism of peak doubling in the superfluid phase. 

The fermion momentum distribution is plotted in Fig. 071. The finite imaginary part of the spectral 
function modifies strongly the momentum distribution. The momentum distributions obtained from the 
spectral functions A and A s are much broader then the one obtained taking the quasi-particle pole approx- 
imation for the spectral function. As expected, the nucleon scattering produces a high momentum tail in 
the distributions. The effect of the superfluid gap can be seen in the difference between the solid and dashed 
lines in Fig. [Tt]. As expected jn], [32| the full spectral function A s gives broader momentum distributions. 
The BCS broadening with small energy gap results in the smearing of the momentum distribution in an 
interval 2A(pf) around the Fermi energy. The density given by the quasi-particle approximation to the 
spectral function A is ~ 1.01 times the actual fermion density given by (p9|). The densities calculated from 
the spectral function A and A s are almost the same. 

In Fig. Hwe plot the temperature dependence of the superfluid gap A(0) as function of the temperature. 
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Figure 17: Nucleoli momentum distribution ob- 
tained from the normal spectral function (solid 
line), from the full spectral function (dashed line), 
and using the quasi-particle pole approximation for 
the spectral function (dotted line) , for the same pa- 
rameters as in Fig. fl2[ 



Figure 18: The superfluid energy gap as function 
of temperature for the BCS theory (dashed line) 
and for the self-consistent calculation (triangles), 
for p — .45(00- 
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Figure 19: BCS quasiparticles energies as function 
of the Hartree-Fock single-particle energy for T = 
1.4 MeV and p = A5p . 



Figure 20: Nucleon momentum distribution for the 
Hartree-Fock single-particles energies (dotted line) 
and the BCS nucleon distribution (solid line), for 
the same parameters as in Fig. nil 
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Below T c = 1.63 MeV the superfluid gap rapidly sets in. We could not perform calculations below T = 1.2 
MeV, due to limited numerical resolution in energy. However, a saturation in the dependence of the 
superfluid gap can already be observed at the lowest temperatures studied. The dashed line denotes the 
result of the BCS theory. The critical temperature for the BCS theory is the same as in the quasi-particle 
approximation, T c = 5.03MeV. The scattering of nucleons and the pseudogap formation above T c make 
the pairing more difficult. In the pairing approximation p4| part of the actual energy gap is due to the 
pseudogap and part to the superfluid gap. Thus the superfluid transition occurs only when the pseudogap 
alone is not enough to exclude the pairing. It occurs at a much lower temperature. Another effect may 
be related to quantitative differences between the quasi-particle and the self-consistent calculations. The 
real part of the self-energy and the T-matrices are different in the two cases. Below T c the self-consistent 
energy gap is always significantly smaller than the BCS result. This is different from the pseudogap pairing 
approximation J24J], where at zero temperature the pseudogap goes to zero and the usual BCS result is 
recovered. We could not perform calculations at low temperatures, due to limitations in computational 
resources. However, already in Fig. 081 a saturation of the energy gap in the self-consistent calculation is 
visible, around A(0) = 7 MeV. 

For comparison we plot also the dispersion relation for the peaks in the BCS spectral function at 



T = 1.4 MeV: E p = ±-i/£| + A 2 (p). Since the energy gap is larger than in the self-consistent calculation 
the splitting of the two energies in the spectral function is much more pronounced. It should be noted that 
the condition of weak BCS, A(0) <C £, P =o, is not really fulfilled in the above case. It is another indication 
that weak coupling BCS theory with free nucleon-nucleon interaction should not be used in the nuclear 
matter. The strong modification of the spectral function by the anomalous self-energy in the BCS case 
is visible also in the fermion momentum distribution (Fig. |20f) . The BCS fermion distribution is much 
broader than the Hartree-Fock distribution. Also the density, which is taken to be p = .45po for the BCS 
fermion distribution is .43po for the Hartree-Fock distribution. 

6 Conclusions 

We have presented a new way of performing nuclear matter calculations in the presence of nuclear superflu- 
idity. The starting point is the T-matrix approximation for the nucleon self-energy. It allows to calculate 
the spectral properties of nucleons in an interacting system, treating on equal footing particle-particle and 
hole- hole diagrams. It leads also to the formation of Cooper pairs below some critical temperature. Strong 
scattering in the vicinity of the critical point, leads to a breakdown of the quasi-particle picture, confirming 
the observation of |rifl . This requires the use of self-consistent, and not quasi-particle, fermion Green's func- 
tions for the calculation of the self-energy A more important generalization of the formalism is required 
to address the superfluid transition. In section such a scheme is presented at least for the case, where 
the T-matrix approach indicates a second order transition. It consists in the use of off-shell propagators in 
the T-matrix ladder, one of them including also the anomalous self-energy. In such a way the superfluid 
gap equation is equivalent to the condition that the T-matrix has a singularity at Fermi energy and zero 
total momentum for all temperatures below T c . The full nucleon spectral function includes the anomalous 
self-energy from the condensed pairs, and the normal self-energy due to nucleon-nucleon scattering. Both 
contributions are important for strongly attractive potentials in nuclear matter leading to an interplay of 
the pseudogap and superfluid gap around T c . The method presented in this work has several important 
properties 



• 



• 



Above T c , it leads to a self-consistent T-matrix resummation, which allows for the description of the 
feedback of the pseudogap formation and of the scattering on the T-matrix. 

It provides a link between the gap equation and the T-matrix equation, at the critical temperature 
and below. 



• It gives the BCS theory in the limit of small scattering rates, unlike the approach Ej. E| 



• 



The value of the energy gap and the critical temperature are strongly reduced, in comparison to the 
BCS (quasi-particle) results. 



When used with realistic nuclear forces, this approach presents a procedure for calculating the nuclear mat- 
ter properties, taking into account self-consistently single-particle spectral properties and fermion pairing. 
Both these important questions were not studied up to now in actual nuclear matter calculations. The first 

17 



of these problems was discussed in a number of works related to high T c superconductivity |23|, gl], p5[ ; 
the first such calculation being presented by Hausmann Qf]. The second question was addressed, to the 
author's knowledge, only by two groups in condensed matter physics |2^, E4J. The work |2|] uses a T-matrix 
resummation with both fermion propagators in the ladder including the anomalous self-energy However, 
due to the truncation of the resummation in the anomalous sector, the usual BCS gap equation is used. In 
that way, the gap equation is not consistent with the T-matrix resummation of the normal self-energy. The 
approach B is similar, i.e. it uses two full propagators in the ladder diagrams. Thus, the gap equation is 
recovered only for small gap, or is used in its linear form. The work p4fl uses the pairing approximation, 
with one mean-field propagator and one with normal and anomalous self-energies. Both self-energies have 
a BCS like form. Our approach is a generalization of this procedure, using two off-shell propagators in the 
T-matrix it allows us to use the off-shell propagator for the calculation of the Hartree-Fock and scattering 
self-energies (O) and (10|) , it is not restricted to a BCS like form of the normal self-energy, and it gives the 
most general T-matrix self-energy scheme above T c . It should be stressed that there is no reason to neglect 
off-shellness in one of the propagators in the T-matrix ladder. The consistency with the gap equation 
requires only to neglect the anomalous self-energy in one of these propagators. The value of the pseudogap 
goes to zero at low temperature in the pairing approximation |24|. On the other, hand the energy gap in 
the self-consistent calculation is much smaller than the BCS gap also at low temperatures (Fig. fiq ). 

The numerical results are given for a very simple interaction. We plan to study the superfluid nuclear 
matter using more realistic separable interactions and also at zero temperature. This would allow to study 
questions like the effect of pairing^ interaction on the ground state energy, the interplay between pairing in 
isospin and 1 channels, and the pairing in neutron matter. Some more fundamental questions remain also 
open. We use for the energy gap a mean-field approximation. Moreover, the interaction potential in the gap 
equation is the free one. Obviously modification of the pairing potential in medium, due to screening, are 
expected. For consistency the same screened potential should be used in the T-matrix equation for normal 
nucleon self-energy, meaning significant complications fl36| . Another, question is related to the cause of the 
reduction of T c . In our calculation both the pseudogap and the nucleon scattering effects are included. 
To disentangle them the self-consistent calculation should be compared to results of the quasi-particle 
approximation and of the pairing approximation below T c . 
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Appendix 

In this appendix we sketch the numerical method used in the iterative solution of the coupled equations for 
the T-matrix, the the self-energy and the spectral function. Special methods must be used for an efficient 
estimation of multidimensional integrals, over energies and momenta, involved in the calculation of the 
T-matrix and the self-energy when using off-shell propagators. 

Let us enumerate below the steps of the iterative procedure, explaining in details the calculation of the 
energy integrals. We present the case of a separable potential of rank 1, as used in the actual calculations 
in this work. In that case the T-matrix in a given channel c is given by 

<k\T°+(P,u)\k'>= ^fJ'fV ( 41 ) 

where the function J c depending only on the total momentum and energy is given by 

A{jp u J -u)')A s (p 2 ,u )(l - f(u> -u>') - f(u )) 



uj — uj + ie 

:>2 



(42) 



Pi,2 = -P74 + fc 2 ± Pk cos(O) 



3 Only Ref. |25[ uses real time formalism for numerical calculations of the self-consistent T-matrix, like the present work. 
4 The effect ofthe off-shellness of the propagators on the ground state energy has not been studied either. 
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For rank one separable interaction, a nonzero solution of the gap equation (|29|) is equivalent to the condition 
1 — A C J C (P = 0,cj = 0) = and the momentum dependence of the pairing gap is A(fc) = A(Q)g c (k)/g c (0). 
Let us present in the following, step by step, one iteration starting from the calculation of the spectral 
function : 

1. Let us assume that the imaginary part of the self-energy is known PL The dispersive contribution to the 
real part of the self-energy can be obtained from (|l4|). The calculation of the spectral functions A and 
A s requires also the knowledge of the chemical potential /i and the pairing gap A. The condition for 
superfluidity is checked and if necessary the gap equation (E9h is solved, using a numerical procedure 
for the solution of a nonlinear equatioifl. For each trial value of A the chemical potential and the 
corresponding Hartree-Fock energy are obtained with the constraint (p9| ) and using Eqs. (|4C|) and 
(|l3|). For each trial value of \x the coupled equations (|13|), J4C|), (g), (|39|) are solved by iteration. As 
a result we obtain the spectral functions A and A s . 

2. In this step we calculate the T-matrix T C (P, ui) for a range of values of total momentum and energy. 
Let us first calculate the imaginary part of the function J c (uj,P) 

l^J c {P,u) = - J ^- J ^ J dcos{Q)A s (p 2 ^'){l- f{u)) 

A{ Pl ,oo - w')(l - /( w - Lo'))g 2 c (k) + ... (43) 

the dots represent a similar term with factors / instead of (1 — /). The above ui integral is a 
convolution integral and can be performed using fast Fourier transform algorithms for numerical 
convolutions. First the two factors in the energy integral are Fourier transformed 

F is) ( P ,t) = FFT[A {s) (p,Lo)] . (44) 

Then the function 

/h 2 rlk f 
— J dco S (Q)F( Pl ,t)F s (p 2 ,t)g 2 c (k) , (45) 

is calculated using standard integration procedures for the two-dimensional momentum integration. 
The first term in Eq. ( |43| ) is obtained by inverse Fourier transform of J. Analogously the second 
term in (|43| ) can be calculated. The real part of J is calculated using a dispersion relation 

„ / n „ f dto —lmJ c (P,u) , „ 

ReJ c (P,w) = V / ^-M. 46 

J IT U) — IjJ 

3. The imaginary part of the self-energy (M) can also be written as a sum of two convolution integrals 
in energy and can be calculated in a similar way as ImJ in the previous point. This gives a next 
iteration of the single-particle width and the iteration returns to the first point. 

A generalization to a separable potential of higher rank is obvious. Since the Fourier transforms are 
performed before the loop over the total momentum, the most expensive numerical cost comes from the 
two dimensional integrals, when calculating J(P, t) for a range of values of t and P. Very similar integrals 
must be performed when calculating the T-matrix and the self-energy in the quasi-particle approximation, 
except that here we have several such integrals. In practice, when starting the iteration from a self-energy 
obtained in a previous calculation at similar temperature and density, around 10 iteration are sufficient to 
converge with a relative deviation ~ 10~ 4 . 

The numerical procedure in its present form introduces a cutoff in momenta of nucleons, to limit the 
energy range. Thus, we cannot treat more realistic potentials, whose cutoff functions g(k) involve large 
momenta. On the other hand, with decreasing temperature, the spectral function becomes narrow near the 
Fermi energy, and cannot be discretized on a finite grid of energy]]. The extension of the present work to 
low energies and to large nucleon momenta requires a combined use of quasi-particle approximation near 
the Fermi energy and for large momenta and a continuum spectral function elsewhere, similar to the ansatz 
usedinRef. 151. 



5 In the first iteration a constant single-particle width or a stored array of ImS obtained earlier for a similar temperature 
and density is used. 

6 This requires the calculation of the inverse T-matrix for a particular value of momentum and energy (the gap equation) 
using a procedure to be described later. 

7 The grid must be equally spaced in order to use fast Fourier transform algorithms 

19 



References 

D.W. Sprung, in Advances in Nuclear Physics, edited by M. Baranger and E. Vogt, (Plenum Press, 
New York, 1972); C. Mahaux and R. Sartor, in Advances in Nuclear Physics, edited by J.W. Ncgele 
and E. Vogt, (Plenum Press, New York, 1991) 

A.L. Fetter and J.D. Walecka, Quantum Theory of Many-Particle Systems (McGraw-Hill, New York, 
1971). 

[3] L.N. Cooper, R.L. Mills and A.M. Sessler, Phys. Rev. 114 (1959) 1377. 

[4] W.H. Dickhoff, Phys. Lett. B210 (1988) 15; B.E. Vonderfecht, C.C. Gearhart, W.H. Dickhoff, A. Polls 
and A. Ramos, Phys. Lett. B253 (1991) 1. 

[5] B.E. Vonderfecht, W.H. Dickhoff, A. Polls and A. Ramos, Nucl. Phys. A555 (1993) 1. 

[6] B.E. Vonderfecht, W.H. Dickhoff, A. Polls and A. Ramos, Phys. Rev. C44 (1991) R1265. 

[7] H.S. Kohler, Phys. Rev. C46 1687 (1992) 1687. 

T. Aim, G. R6pke, A. Schncll, N.H. Kwong and S. Kohler, Phys. Rev. C53 (1996) 2181; A. Schncll, 
T. Aim and G. Ropke, Phys. Lett. B387 (1996) 443. 

[9] D.J. Thouless, Ann. Phys. (NY.) 10 (1960) 553. 

[10] L.P. Kadanoff and P.C. Martin, Phys. Rev. 124 (1961) 670. 



[11] P. Bozek, aucl-th/9811073 



[12] F. de Jong and H. Lenske, Phys. Rev. C56 (1997) 154. 

[13] M.I. Haftel and F. Tabakin, Nucl. Phys. A158 (1970) 1; M. Golberger and K. Watson, Collision 
Theory (John Wiley and Sons, New York, 1964). 

[14] L.P. Kadanoff and G. Baym, Quantum Statistical Mechanics (Benjamin, New York, 1962). 

[15] P. Danielewicz, Ann. Phys. 152 (1984) 239. 

[16] W. Botcrmans and R. Malfliet, Phys. Rep. 198 (1990) 115. 

[17] M. Schmidt, G. Ropke and H. Schulz , Ann. Phys. 202 (1990) 57. 

[18] Y. Yamaguchi, Phys. Rev. 95 (1954) 1628. 

[19] M. Baldo, I. Bombaci, G. Giansiracusa, U. Lombardo, C. Mahaux and R. Sator, Nucl. Phys. A545 
(1992) 741. 

[20] O. Benhar, A. Fabrocini and S. Fantoni, Nucl. Phys. A550 (1992) 201. 

[21] M. Randeria, N. Trivedi, A. Moreo and T. Scalettar, Phys. Rev. Lett. 69 (1992) 2001; P.G. McQueen, 
D.W. Hess and J.W Serene, Phys. Rev. B 50 (1994) 7304; N. Trivedi and M. Randeira, Phys. Rev. 
Lett. 75 (1995) 312; R. Micnas, M.H. Pedersen, S. Schafroth, T. Schneider, J.J. Rodriguez-Ninez and 
H. Beck, Phys. Rev. B52 (1995) 16223; MY. Kagan, R. Fresard, M. Capezzali and H. Beck, Phys. 



Rev. B57 (1998) 5995; M. Randeira, cond-mat/9710223 



[22] M.H. Pedersen, J.J. Rodriguez-Ninez, H. Beck, T. Schneider and S. Schafroth, Zeit. fur Phys. B103 
(1997) 21. 



[23] J. Maly, B. Janko and K. Levin, cond- mat/9805018; J. Maly, B. Janko and K. Levin Phys. Rev. B56 



(1997) R11407; I. Kosztin, Q. Chen, B. Janko and K. Levin, Phys. Rev. B58 (1998) R5936. 
[24] Q. Chen, I. Kosztin, B. Janko and K. Levin, Phys. Rev. Lett. 81 (1998) 4708. 
[25] B. Kyung, E. G. Klepfish and P.E. Kornilovitch, Phys. Rev. Lett. 80 (1998) 3109. 
[26] R. Hausmann, Phys. Rev. B49 (1994) 12975. 

20 



[27] R. Hausmann, Zeit. Fiir Phys. B91 (1993) 291. 

[28] A.M. Sessler, in Liquid Helium, Proceedings of the International School of Physics Enrico Fermi, edited 
by G. Careri, (Academic Press, New York, 1963); V.G. Valeev, G.F. Zharkov and Yu.A. Kukharcnko, in 
Nonequilibrium Superconductivity, Proceedings of the Lcbcdcv Physics Institute vol. 174, edited by V.L. 
Ginzburg, (Nova Science Publishers, New York, 1988); D. Rainer and J. A. Sauls, in Superconductivity, 
Lecture Notes of the ICTP Spring College in Condensed Matter on "Superconductivity" , edited by 
P.N. Butcher and Yu Lu, (World Scientific, Singapore, 1995). 

[29] A.L. Goodman, Nucl. Phys. A186 (1972) 475. 

[30] M. Baldo, U. Lombardo and P. Schuck, Phys. Rev. C52 (1995) 975. 

[31] A.B. Migdal, Theory of Finite Fermi Systems (John Wiley and Sons, New York, 1967). 

[32] J.R. Schricffer, Theory of Superconductivity (W.A. Benjamin, Inc., Massachusetts, 1964); G.D. Mahan, 
Many-Particle Physics (Plenum Press, New York, 1981). 

[33] P.W. Anderson and P. Morel, Phys. Rev. 123 (1960) 1911. 

[34] R. Balian and N.R. Werthamcr, Phys. Rev. 131 (1963) 1553. 

[35] P.W. Anderson and W.F. Brinkman, in The Physics of Liquid and Solid Helium, part II, edited by 
K.H. Benneman and J.B. Ketterson, (John Wiley and Sons, New York, 1978). 

[36] A.D. Jackson, A. Lande and R.A. Smith, Phys. Rep. 86 (1982) 55. 



21 



